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On the Synchronization of Second-Order Nonlinear Systems 
with Communication Constraints 
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Abstract 

This paper studies the synchronization problem of second-order nonlinear multi-agent sys¬ 
tems with intermittent communication in the presence of irregular communication delays and 
possible information loss. The control objective is to steer all systems’ positions to a common 
position with a prescribed desired velocity available to only some leaders. Based on the small- 
gain framework, we propose a synchronization scheme relying on an intermittent information 
exchange protocol in the presence of time delays and possible packet dropout. We show that 
our control objectives are achieved with a simple selection of the control gains provided that the 
directed graph, describing the interconnection between all systems (or agents), contains a span¬ 
ning tree. The example of Euler-Lagrange systems is considered to illustrate the application 
and effectiveness of the proposed approach. 


1 Introduction 

Motion coordination of nonlinear multi-agent systems has received an increased interest in the 
control community due to the potential applications involving groups of robotic systems and au¬ 
tonomous vehicles in general m ■ Multi-agent systems control can be formulated as synchro¬ 
nization or consensus problems, where the goal is to drive the networked systems (or agents) to a 
common state using local information exchange. Other related problems include flocking, swarming, 
and formation control of mechanical systems. Built around the solutions of the consensus problem 
of linear multi-agent systems, several coordinated control schemes have been recently developed for 
second-order nonlinear dynamics, which can describe various mechanical systems, with a particular 
interest to leaderless synchronization problems E®, cooperative tracking with full access to the 
reference trajectory EH®, leader-follower with single leader &m or multiple leaders EMSI, 
to name only a few. Algebraic graph theory, matrix theory, and Lyapunov direct method have 
been shown useful tools to address various problems related to the systems dynamics, such as 
uncertainties, and the interconnection topology between the team members. 

In addition, various recent papers address the synchronization problem of nonlinear systems by 
taking into account delays in the information transfer between agents, which is generally performed 
using communication channels. In [7] and m, it has been shown that output synchronization 
of nonlinear passive systems is robust to constant communication delays if the interconnection 
graph is directed, balanced and strongly connected. A similar property was shown in 0 under 
unbalanced directed graphs using the contraction theorem. In [15] , a delay-robust control scheme 
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is proposed for relative-degree two nonlinear systems with nonlinear interconnections. With the 
same assumption on the delays, adaptive synchronization schemes have been proposed in )19H20| 
for networked robotic systems under a directed graph. In addition to constant delays, a virtual 
systems approach has been suggested in pT][22] to account for input saturations and to remove the 
requirements of velocity measurements. Control schemes that consider time-varying communication 
delays have also been proposed for the attitude synchronization of rigid body systems [251121]. 
formation control of unmanned aerial vehicles |25| , and consensus of networked Lagrangian systems 
]26] . yet in the case of undirected interconnection graphs. More recently, a small-gain framework is 
proposed in m for the synchronization of a class of second-order nonlinear systems in the presence 
of unknown irregular time-varying communication delays under general directed interconnection 
topologies. 

One important problem when dealing with second-order nonlinear systems in the presence of 
communication delays is to achieve position synchronization, i. e ., all positions converge to a common 
value, with some non-zero final velocity. In fact, in most of the above mentioned synchronization 
laws with communication delays, a static leader or no leader are assumed and position synchroniza¬ 
tion is achieved with zero final velocity. The only cases where the final velocities match a non-zero 
value assume a full access to a reference trajectory or to a leader’s states (position and velocity). 
By full access, it is meant that this information is available to all agents without delays. The main 
challenge in this case resides in the fact that imposing a non-zero final velocity ultimately requires 
some information on the delays to achieve position synchronization. In fact, a possible solution to 
this problem might be to explicitly incorporate the delays in the control algorithms as suggested 
in 1281129] for linear second-order multi-agents. This, however, comes with the assumptions of full 
access to the desired velocity and the communication delays are exactly known. 

Another issue that can be observed in all the aforementioned results is the assumption that 
information is transmitted continuously between agents. In fact, it is not clear if these results 
still apply in situations where agents are allowed to communicate with their neighbors only during 
some disconnected intervals (or at some instants) of time. This can be induced by environmen¬ 
tal constraints, such as communication obstacles, temporary sensor/communication-link failure, 
or imposed to the communication process to save energy/communication costs in mobile agents. 
For linear first-order multi-agent systems, the authors in [30] have proposed a consensus algorithm 
based on the output of a zero-order-hold system, which is updated at instants when the informa¬ 
tion is received and admits as input the relative positions of interacting agents. In the presence of 
sufficiently small constant communication delays and bounded packet dropout, the proposed dis¬ 
continuous algorithm in m achieves consensus provided that self-delays are implemented and the 
non-zero update period of the zero-order-hold system is small. A similar approach has been applied 
for double integrators in [31 [l32l . where asynchronous and synchronous updates of the zero-order- 
hold systems have been addressed, respectively, without communication delays. Here, synchrony 
means that all agents receive information at the same instants. In [53], a switching algorithm has 
been proposed for second-order multi-agents in cases where communication between agents is lost 
during small intervals of time, yet without communication delays. The latter result has been ex¬ 
tended to multi-agent systems with general linear dynamics [34] and globally Lipschitz nonlinear 
dynamics [35] , where it has been shown that consensus can be achieved under some conditions on 
the communication rates and interaction topology. 

In this paper, we consider the synchronization problem of a class of second-order nonlinear 
systems with intermittent communication in the presence of communication delays and possible 
packet loss. Here, it is required that all systems achieve position synchronization with some non¬ 
zero desired velocity available to only some systems in the group acting as leaders. Based on the 
small-gain approach, we propose a distributed control algorithm that allows agents to communicate 
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with their neighbors only at some irregular discrete time-intervals and achieve our control objec¬ 
tive. A discrete-time consensus algorithm is also used to handle the partial access to the desired 
velocity. In the case where no desired velocity is assigned to the team, the proposed synchroniza¬ 
tion algorithm achieves position synchronization with some velocity agreed upon by all agents. In 
both cases, it is proved that, under some sufficient conditions, synchronization is achieved in the 
presence of unknown irregular communication delays and packet loss provided that the intercon¬ 
nection topology between agents is described by a directed graph that contains a spanning tree. 
The derived conditions impose a maximum allowable interval of time during which a particular 
agent does not receive information from some or all of its neighbors. This interval, however, can be 
specified arbitrarily with a choice of the control gains. To illustrate the applicability of the proposed 
approach, we derive a solution to the above problems in the case of networked Lagrangian systems, 
and simulation results that show the effectiveness of the proposed approach are given. 

2 Background 

2.1 Graph theory 

Let G = (A f,£) be a directed graph, with a set of nodes (or vertices) Af , and a set of ordered 
edges (pairs of nodes) £ C Af x Af. An edge (j, i) £ £ is represented by a directed link (arc) 
leaving node j and directed toward node i. A directed graph Q is said to contain a spanning tree if 
there exists at least one node that has a “directed path” to all the other nodes in the graph; by a 
directed path (of length q) from j to i is meant a sequence of edges in a directed graph of the form 
(j, Zi), , (lq-i,lq), with l q = i, where for q > 1 the nodes j , Zi ,..., l q - 1 £ Af are distinct. 

Node r is called a root of Q if it is the root of a directed spanning tree of G\ in this case, Q is said 
to be rooted at r. 

Given two graphs Gi = (AT, £i), G2 = (AT, £2) with the same vertex set AT, their composition 
G3 '■= {Af, £3) = Gi o G2 is the graph with the same vertex set AT where (j, i) £ £3 if and only if 
( j , l) £ £2 and (l,i) £ £\ for some l £ Af. Composition of any finite number of graphs is defined by 
induction. In the case where Gi and G2 contain self-links at all nodes, the edges of Gi and G2 are 
also edges of Q3. In this case, the definition above also implies that G3 contains a path from j to i 
if and only if G2 contains a path from j to l and Gi contains a path from l to i. A finite sequence 
of directed graphs Gi , G2, ..., Q q with the same vertex set is jointly rooted if the composition 
G q 0 Gq-i o ... o Gi is rooted. An infinite sequence of graphs Go , G 1, ... is said to be repeatedly 
jointly rooted if there exists k* £ such that for any a £ Z + the finite sequence Go , G0+1, ■ ■ 
Go+k* is jointly rooted. (See [ 3 ja] for more details on graph composition). 

A weighted directed graph Q w consists of the triplet (Af , £, A), where Af and £ are, respectively, 
the sets of nodes and edges defined as above, and A is the weighted adjacency matrix defined such 
that an = 0, a-y > 0 if ( j,i ) £ £, and aij = 0 if (j, i) ^ £. Note that thus defined graph does not 
contain self-links at any node and will have the same properties as the unweighted graph with the 
same sets of nodes and edges. The Laplacian matrix L := [Uj] £ M nxn of the weighted directed 
graph Gw is defined such that: la = )Cj=i a ij> an( i hj = ~ a ij f° r * 7 ^ j- 

2.2 Stability Notions and preliminary result 

Consider an affine nonlinear system of the form 

x = f(x) + Y^ =1 9i(x)ui, 

Vi = hj{x), j = l,...,q, 


3 


where x £ M w , ut £ R mi for i £ Af p := {1,... ,p}, yj £ R m > for j £ Af q := {1,... ,q}, and /(•), gi{-), 
for i £ Afp , and hj(-), for j £ Af q , are locally Lipschitz functions of the corresponding dimensions, 
/(0) = 0, h( 0) = 0. We assume that for any initial condition x(to) and any inputs u\(t), ..., u p (t) 
that are uniformly essentially bounded on \ta,t\), the corresponding solution x(t) is well defined 
for all t £ [to, h]. 

Definition 1 . JS'A/ A system of the form (JT]) is said to be input-to-state stable (ISS) if there exisf] 
(3 £ /Coo o-nd 7 j £ 1C, j £ Af p , such that the following inequalities hold along the trajectories of the 
system for any Lebesgue measurable uniformly essentially bounded inputs Uj, j £ Af p : 

i) uniform boundedness: V to, t £ K., t > to, we have 

\x(t)\ < /3(|a;(t 0 )|) + E Ij ( sup M«)l ) , 
j =l / 


ii) asymptotic gain: 

p / 

limsup |ar(t)| < E Ij ( limsup |uj(t)| 

t —^ -\- OO j=l y t —^ oo 

In the above definition, | • | denotes the Euclidean norm of a vector and 7 j £ 1C, j £ Af p , are 
called the ISS gains. It should be pointed out that for a system of the form ©, the ISS implies the 
input-to-output stability (IOS) [37], which means that there exist Pi £ ICC and 7 y £ 1C, i £ Af q , 
j £ Afp , such that the inequality 


< Pi (|a;(fo)|,f) + E 7 ij ( sup |uj(s)| ] , 
j— 1 \se[to,t) / 

holds for all i £ Af q and V to, t £ M, t> to. In this case, the function 7y £ 1C, i £ Af q and j £ Af p , is 
called the IOS gain from the input Uj to the output yi. In the subsequent analysis, we will mostly 
deal with the case where the IOS gains are linear functions of the form "fij(s) := 7 ?- • s, where 
lij > 0; in this case, we will simply say that the system has linear IOS gains 7 ?- > 0. 

The convergence analysis in this paper is based on the following small-gain theorem. 

Theorem 1. Consider a system of the form ©. Suppose the system is IOS with linear IOS gains 
1% > 0. Suppose also that each input Uj(-), j £ Af p , is a Lebesgue measurable function satisfying: 
Uj(t) = 0, for t < 0, and 


\uj(t)\ < £ N i ■ sup |y»(s)| + \Sj(t)\, ( 2 ) 

i e jy q 

for almost all t > 0, where fiji > 0, all dji(t ) are Lebesgue measurable uniformly bounded nonnega¬ 
tive functions of time, andSj(t) is an uniformly essentially bounded signal. LetT := T°• AA £ R 9XlJ , 
where T° := {7 ?•}, AA := {pji}, i £ Af q , j £ Af p . If p(T) < 1, where p(r) is the spectral radius of 
the matrix T, then the trajectories of the system m are well defined for all t > 0 and such that all 
the outputs gift), i £ Af q , and all the inputs Uj(-), j £ Af p , are uniformly bounded. If, in addition, 
|<5j(t)| —>■ 0 at t —> +oo ; j £ Afp, then \yi(t)\ —> 0, \uj(t)\ —>0 as t —► +00 for i £ Af q and j £ Af p . □ 

Theorem [T| is a version of EH Theorem 1 ] and is also a special case of the result given in [3S]; in 
particular, its proof follows the same lines as in the proof in [271 Theorem 1], and hence, is omitted. 

1 The definition of the class functions /C, /Coo, and ICC can be found in |38| . Also, 1C := /CU{0}, 1C oo := /CooU{0}, 
where O is zero function, O(s) = 0 for all s > 0. 
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3 Problem Formulation 


3.1 System Model 


Consider n not necessarily identical second-order nonlinear systems (or agents) governed by the 
following dynamics 


Pi{t) = Vi(t), 

Vi(t) = Fi(pi(t),Vi(t),Ui(t)), 


(3) 


where pi £ R m and v\ are the position-like and velocity-like states, respectively, tq £ R m are the 
inputs, and A f := {1,..., n}. The functions Fj are assumed to be locally Lipschitz with respect to 
their arguments. Note that equations <01 can be used to describe the full or part of the dynamics 
of various physical systems. 

The systems m are interconnected in the sense that some information can be transmitted 
between agents using communication channels. This interconnection is represented by a directed 
graph Q = (A/",£), where A f is the set of all agents, and an edge (j,i) £ £ indicates that the i-th 
agent can receive information from the j-th agent; in this case, we say that j and i are neighbors 
(even though the link between them is directed). While the interconnection graph Q is fixed, 
the information exchange between agents is not continuous but discrete in time and is subject to 
communication constraints as described in the next subsection. 


3.2 Communication Process 

In this paper, we consider the case where the communication between agents is intermittent and is 
subject to time-varying communication delays, information losses, and blackout intervals. Specifi¬ 
cally, it is assumed that there exists a strictly increasing and unbounded sequence of time instants 
tk '■= kT £ R_|_, k £ Z + = {0,1,...}, where T > 0 is a fixed sampling period common for all agents, 
such that each agent is allowed to send its information to all or some of its neighbors at instants tk, 
k £ Z + . In addition, for each pair (j, i) £ £, suppose that there exist a sequence of communication 
delays (r^’^)fcgz + that take values in {R + U +oo} such that the information sent by agent j at 
instant tk can be available to agent i starting from the instant tk -)-tP’^ . In particular, it is possible 
that = Too for some k £ Z + , which corresponds to a situation where agent j has not sent 

information at instant tk to neighbor i at all, or the corresponding information was never received 
possibly due to packet loss in the communication channel. The following assumption is imposed on 
the communication process between neighboring agents. 

Assumption 1. For each ( j,i ) £ £, there exist numbers k* £ N, h > 0, and an infinite strictly 
increasing sequence := {kQ ,l \ kff’ l \ ...} C {0,1,...} satisfying 

i) kg' 1 ^ < k*, and k\^ — < k*, l £ {0,1,.. 

ii) < h for each k £ 

Assumption Q] essentially means that, for each pair (j,i) £ £, and per any k* consecutive 
sampling instants, there exists at least one sampling instant at which agent j has sent information 
to agent i, and this information has been successfully delivered with delay less than or equal to h. 
Note that h is not an imposed upper bound of the delays; in particular, the possible case where 
Ta’^ > h for some a ^ Kfj,*) is not excluded. Assumption 0 also implies that, for each pair 
( j , *) £ E , the maximal interval between two consecutive instants when agent i receives information 
from agent j is less than or equal to 

h*:=k*T + h. (4) 
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It is worth pointing out that k*, h , and h* are considered common to all (j,i) £ £ for simplicity, 
and can be seen as the maximum of the corresponding parameters defined for each pair (j,i) £ £. 
Also, the above assumption does not require that all agents broadcast their information to some or 
all of their prescribed neighbors at the same instants of time. 

3.3 Problem statement 

Consider multi-agents © interconnected according to Q and the communication between agents 
satisfies Assumption |I] Suppose also that a constant desired velocity Vd £ R m is available for a 
subset C C AT of agents, called leaders. The rest of the systems (belonging to the complementary 
subset T := Af\C) are referred to as followers. Our goal is to design synchronization schemes for 
the nonlinear multi-agent system © such that the following objectives are attained. 

Objective 1. In the case C ^ 0, it is required that Pi(t) — Pj(t) —> 0, Vi(t) —> Vd, as t —»• -foo for 
all i,j £ A I. 

Objective 2. In the case C = 0, it is required that Pi(t) — Pj(t ) —► 0 and Vi(t) —> v c as t —► +oo 
for all i,j £ A f and for some final velocity v c £ R m . 

To achieve the above objectives, we adopt an approach that takes its roots from the control 
of robotic systems |3D] and has been recently used to address various synchronization problems of 
mechanical systems (see, for instance, [5l HTllT9l[25li27] I . To explain this approach, let v ri (t) £ C 1 be 
a reference velocity for the i-th system, for i £ A f. Equations © can be rewritten in the following 
form, 


Pi(t) = v r i(t) + ei(t), (5) 

ei(t) = Fi(pi,ei + v ri ,uf) - v ri , (6) 

i £ A/”, where := Vi — v Ti is the velocity tracking error. Equations © describe the dynamics of 
agents in a multi-agent system with v Ti (t) being the reference input, while ej(i) is a perturbation 
term with dynamics described by ©■ The synchronization problem can now be solved using a two 
stages approach described as follows. In the first stage, the input u;(i) in © is designed to guarantee 
the convergence of the error signals e»(i) to zero. In the second stage, appropriate algorithms for 
v ri (t) are designed using the position-like states of the systems such that the trajectories of the 
dynamic systems © satisfy Objectives [I] or [2] As mentioned in the Introduction, the problem of 
designing such algorithms for v Vi (t) in the presence of the communication constraints described in 
Section [3.21 is yet unsolved, even in the case where no perturbation term exists, i.e., = 0. In the 

presence of nonzero signals (t), the problem becomes more complicated since there typically exists 
some coupling between the signals ei(t) and v n (t) as they both depend on the states of ©-©. 
Note that, the first step, the design of the control law m in © that guarantees desirable properties 
of the error signal e*, for a given v ri and v ri , can be achieved using various existing approaches to 
tracking control design for nonlinear systems. This step is not addressed in this work; instead, the 
following assumption is made. 

Assumption 2. For each system in © and a given reference velocity signals v ri (t) £ C 1 and 

v ri {t), there exists a static or dynamic tracking control law Ui such that the following hold: 

• the error signal ei{t) is uniformly bounded; 

• if y ri(t) and v ri (t) are globally uniformly bounded, then ei(t) — > 0 as t —» +oo. 

Under Assumption [2] the synchronization problem of the nonlinear multi-agent system © is 
reduced to the design of the reference velocities v ri £ C 1 , * £ Af, such that v ri is well defined 
(available for feedback) and the trajectories of © satisfy Objectives [T| or [2] 
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4 Synchronization Scheme Design 

In this section, we present a method for the design of the reference velocities v ri (t), i £ Af, such 
that Objectives |T] and [2] are achieved using intermittent communication between the agents in the 
presence of time-varying communication delays and information losses. For this, we let p*■(&), for 
each (j,i) £ £ and each k £ Z+, denote the information that can be transmitted from agent j to 
agent i at t = kT. Specifically, p) (k) := [pj ( kT ), Vd, ( k ), k] , where pj ( kT ) is the position-like state of 
the j-th system at t = kT , v ( i :j (k) is a discrete-timqj estimate of the desired velocity obtained by the 
j- th agent according to an algorithm described below, and k is the time-stamp, i.e., the sequence 
number at which information was sent. Also, for each pair (j, *) £ £ and each time instant t > 0, 
we let kij(t ) be the largest integer number such that p*(fc.y(f)) = \pj(kij(t)T), Vd^ (kij(t)), kij(t)] is 
the most recent information of agent j that is already delivered to agent i at t , i.e., 

kij(t) := max{fc £ Z + : kT + rjf’^ < t}. (7) 

Note that the number kij(t) can be determined by a simple comparison of the received time stamps. 
Now, for each i £ Af , the reference velocity v ri ( t ) is designed in the following form, 


v ri{t) = Vi(i) + Vdi(t), (8) 

where Vd t is a sufficiently smooth estimate of the desired velocity available for the i -th agent, and 
rji is a synchronization term designed with the purpose of position synchronization between agents. 
The design of Vdi and 77 i are addressed below in detail. 

4.1 Desired velocity estimation 

In this subsection, we present a method for the design of Vd, hr (JHJ. As explained above, each leader 
has direct access to the desired velocity Vd £ R m . The followers, on the other hand, do not have 
direct access to the desired velocity Vd\ instead, they estimate it through the following discrete-time 
consensus algorithm that is updated at instants aT , for <7 £ Z+, 

( a ) = v d for i £ C, (9) 

VdM + l) = 1 for is J, ( 10 ) 

I T a ) I jeNi(a) 


where 


Vdi } (cr) 


Vdj(kij(crT)) if j^i, 
Vdi(cr) if j = i. 


In the above algorithm (fTUl) - (fTTl) . Ni(a) := {*} U Nij(a ), with 


( 11 ) 


Nij(a) := {j : (j,i) £ £, k i:j {aT) > %((<7 - 1)T)}, 

and | Ni(a) | denotes the number of elements in Ni(cr). Recall from 0 that the vector 7 )d j ( kij (crT)) 
is the most recent desired velocity estimate (obtained by agent j) that is already available to agent 
i at instant aT. As such, the set Nij(a) denotes the set of the neighbors of the i-th follower such 
that the most recent data from these neighbors has been received during the interval ((a — 1)T, aT]. 

2 For simplicity, we use throughout the paper the notation x(k). k G Z+, instead of x(kT) for the discrete-time 
signals. 
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Therefore, the update law for the followers (HD-(EH) is based on each agent’s own estimate of the 
desired velocity as well as on the most recent velocity estimates (received from its neighbors) that 
have not been already used in the update law at earlier sampling instants. The leaders, on the 
other hand, do not update their estimates as can be seen from ©. It should be noted that the set 
Ni(a) used in EHD-EH can be obtained by simple comparison between the successfully received 
time stamps. 

Using the discrete-time desired velocity estimates obtained by each agent, we let 

Vdi = Vd for i £ C , ( 12 ) 

Vdi = —kfidi ~ fcf (vd. ~ Vdi for * £ IF, (13) 

where |_aj denotes the integer part of a € K+, and k\, kf are strictly positive scalar gains. Note 
that Vdi £ C 1 . 


Proposition 1. Consider the discrete-time consensus algorithm (j9]) - (llll) . Suppose Assumption Q] 
holds. If C 7 ^ 0 and the directed interconnection graph Q is rooted at r £ C, then i>di(cr) —»• v c = Vd 
as a —> +00 for all i £ Af. Also, if C = 0 and Q is rooted, then id, (<r) —> v c as a —»• +oo for all 
i £ Af, for some v c £ R m . In both cases, Vdi(t) and Vdi(t ) are uniformly bounded and Vdi(t) —*■ v c 
as t —>■ + 00 , for all i £ AI. □ 


4.2 Design of the synchronization terms 

In this subsection, the design of the synchronization term Vi in (© is addressed. For this, let 
A = [a»j] £ TZ nxn be an arbitrary weighted adjacency matrix (Defined in Section 12.11) assigned 
to the graph Q; the resulting weighted directed graph is denoted by Q w = (Af,£,A). Let Aft := 
{j £ Af : ( j,i ) £ £} denote the set of neighbors of agent i in Q w . Also, let ft,; := V : u ;/ and 
Af i := {i £ Af: ft,^0} denote the subset of all agents that have at least one incoming link in Q w . 
Consider the following design of the synchronization term rji in ((5J) 


Vi = -kiVi- Ai(pj-^i), 

'•pi = -A + V di + j- E jeMi "ul'.i" (*)> 


(14) 


for i £ Af^ 


Vi = Vi = 0 for i £ Af \ Af ^, 


(15) 


where kf , \ are strictly positive scalar gains, Vd, is defined in E2D-E3D, and the vector (t ), for 
all j £ Af i: is an estimate of the current position of the j-th agent defined using the most recent 
information available to the i-th agent at t as 


Pj\t) := Pj(hj(t)T) + eij(t), (16) 

e»jW : = v dj (kij(t)) ■ (t — kij(t)T). (17) 


with kij(t) being defined in Q). 

Remark 1 . Note that, due to the intermittent and delayed nature of the communication process, we 
have considered a dynamic design for the synchronization terms vi- This guarantees that Vi £ C 1 , 
which is difficult to realize using static synchronization terms in view of the irregularities of the 
information received by each agent. In addition, the vectors v i are designed as in E3MH0) such 
that the closed loop system © with © and El-EH) for each agent is IOS with arbitrary IOS gains. 


As will be made clear in the next subsection, by employing the small gain theorem, Theorem[ 7J we 
will show that our control objectives are achieved under some conditions that can be always satisfied. 

4.3 Main result 

The following theorem describes the conditions under which Objectives [T| and [5] are achieved. 

Theorem 2. Consider the network of n systems 0, where the interconnection topology is de¬ 
scribed by a directed graph Q = (A/ - , £) and the communication process between the systems satisfies 
Assumption ]! J Suppose that each system is controlled by a control law Ui satisfying Assumption [1 
where the corresponding reference velocity v ri (t) is generated by (j8j), with (l9l)- (TT3l) and lfl4l) - (fT7l) . 
Let the control gains be selected such that 

Pi>l + 2-h*, forieAf (18) 

where h* > 0 is defined by ©, and 

Pi := - max {JLe{p it i), Tle{pi ^)), 

where pip, pi^ are the roots of x 2 + kfx + Xi = 0. Then, for arbitrary initial conditions, we have 

• Objective Q] is achieved if Q contains a spanning tree with a root r £ C. 

• Objective [H is achieved if Q contains a spanning tree. 

Proof. First, it should be noted from (0, (I9l)- (fl3l) . and (fl4l) - (fl7l) that v Ti £ C 1 , and v ri (t) = rji+Vd, 
can be obtained from the solution of the dynamic systems m and m, and is available for 
feedback. Therefore, applying the control law that satisfies Assumption [2] in 0 guarantees that 
the velocity tracking error efit) = vfit) — v Ti {t) is uniformly bounded. 

For each i £ A, let 


Pi := Pi - ifi, ( 19 ) 

Using the relation e* = (vi — r/i — Vd J with 01 , one can write 

Vi = - Kvi - A ifii 
Pi = Vi +A + o 
ifi = Si T (fi (22) 

for i £ A^, with 



£i = ei + — a n(Pj ~pf)- ( 24 ) 

Ki ieM 


( 20 ) 

( 21 ) 


We can verify that each system m-m with output Vi is IOS with respect to the input vectors 
<f>i and £j. This follows by noticing that the following estimates 


4>i(t) 






+ sup 


kite) I 


sup l&OOl' 
se[t 0 ,t] 


(25) 
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(26) 


Vi{t) 

Pi(t ) 

+ — SUp |£i(?)|, 

Pi «e[t 0 ,t] 


< e 




Vi(t o) 
Pz(io) 


1 


H-sup 

fa(fa 

Pi «6[to,t] 



hold for all t > to, where /q is defined in Theorem [2j More precisely, inequality (l25l) indicates that 
system E121) is ISS with respect to the inputs £.; and fa, with unity ISS gains. Also, ESI) implies 
that EIB - ETl) is ISS with respect to the inputs fa and fa, with ISS gains both equal to Since 
the cascade connection between two ISS systems is ISS, we conclude that E0l) -(|22 |) is ISS. As a 
result, the system m-m with output r]i is IOS with respect to the inputs £j and 4 >i with linear 
IOS gains equal to |. and A-, respectively. 

Therefore, all the systems (20-©, for i G Af^, can be regarded as a system with n ** outputs, 
given by yi = rji, i G Aft, and 2 rfi inputs that can be ordered as: u-n := e%, U2i-i '■= fa, for 
i G Af^ , where nf := | Af ** | denotes the number of elements in A ft. From the above analysis, we 
can conclude that such system is IOS, with IOS gain matrix T 0 := { 7 ° } £ u ™ 8 * 2 " 11 given by 


[ 1 1 Pi 
1u = { 2 //x* 


if l = 2i — 1, i G A 
if l = 2i, i G Aft, 


0 otherwise. 


Moreover, using E3ll - E4j) with (fl6l) - (fl5l> and the relation Vi = (a + rji + v di ) for i G Af, one can 
write 


l«*-i(*)l < E — k'WI + l^-rWI, (27) 

jeM} 

\u 2 i(t)\ < |ej(t)| + E — I \vfas)-v dj {kij(t))\ds 

< E Cll ° h sup \Vj(fa\ + \^2i(t)\, (28) 


where we used Assumption |T] and Q to conclude that (t — kij{t)T) < ( k*T + h) := h* , the set 
Aff := {j G A ft : (j, i) G £} is used here due to csi), and 


\fai-i{t)\=\ei(t)\+ E — | ej(t) + v dj (t) - v di (t)\ 

ieM Ki 


(29) 


M*)l = e 


an ■ h* 


sup | v dj (fa - V dj (kij (f)) | 


K i «e[ky(t)r,t] 

+ \ e i(fa\ + E a, \. h SUP l e J'(OI ■ 

jeA f. se[fcy(t)T,i] 


(30) 


Therefore, one can conclude that the input vectors Uj, j G {1,..., 2 nf}, satisfy the conditions of 
Theorem [TJ where the elements of the interconnection matrix A4 := {yij} G 


rdxrd are obtained 
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as 



otherwise, 


ifi = 2f-l, j&Nl i6 UK 
if l = 2 i, j G A/f, i G A/ 1 , 


and where S-ii-iit) and 821 (t), for i G A ft, satisfy (EHl) - (l30l) . Note that in view of Assumption [5] and 
the result of Proposition [TJ we have 8 j(t), j = {1,..., 2n**} are uniformly bounded. 

Therefore, the elements of the closed-loop gain matrix T := T 0 • A4 = {o'iy in Theorem Q] 

can be written as 



Taking into account the fact that yu = 0 and the elements of T are nonnegative, one can conclude 
using Gersgorin disk Theorem that p(T) < 1 if 7 ij < 1, for i G AA Noting that 7 ij = 

SjeN® jT7F (-*- + 2 • h*) and Af- C A/i, the condition p(F) < 1 is satisfied by (1T5|) . 

Therefore, all the conditions of Theorem [T] are satisfied and one can conclude that r]i, <j>i and 
Si, for i G A - *, are uniformly bounded. In addition, the ISS property of (l2(l - (l22l) guarantees that 
r ].- L , pi and ipi, for i G AT\ are uniformly bounded. This with the result of Proposition [Tj and (fl5l) 
lead to the conclusion that v Ti and v ri , i G A f, are uniformly bounded, and hence Assumption [2] 
guarantees that d{t ) —> 0 as t —> + 00 , for i G Af. Furthermore, using the result of Proposition |T| 
and the fact that kij(t) —> +00 as t —> + 00 , it can be verified from (l^l) - (l301) that: Sj(t) —> 0 at 
t + 00 , j G {1,..., 271**} if Q is rooted at r G C ^ 0, or Q is rooted in the case £ = 0. 

Consequently, one can conclude from Theorem [T| that rji(t) 0, </>i(t) —> 0, £j (t) —> 0 for i G A 
as t —> + 00 . Since (rji(t) + ei(t)) = (Vi(t ) — v^it)) —> 0 as t —> + 00 , for i G AT, the result of 
Proposition [T] implies that Vi(t) —> Vd as t —> +00 if Q is rooted at r G £ ^ 0. The same proposition 
implies that Vi(t) —> v c as t —> + 00 , for some v c G if Q is rooted and £ = 0. 

In addition, since system ([20]) - (f22|l is ISS, we have pi(t) —>■ 0 and —> 0 as t —> +00 

for i G A/”*. Using (TTOli . one gets Ki{pi + ipi) = X^j=i a b(£* — Pj) i s uniformly bounded and 
1 a ij(Pi(t) ~ Pj(t)) —> 0 as t —> +00 for all i G A fK This with the fact that m = 0 for i G Af\AT$ 
lead to the conclusion that (L ® l m )p(t) —> 0 as t —> + 00 , where L is the Laplacian matrix of the 
interconnection graph Q w , I m is the m x m identity matrix, p G R” m is the vector containing all pt 
for i G Af , and ® is the Kronecker product. Finally, since (L ® I m )p = 0 implies that p\ = ... = p n 
if Q w contains a spanning tree 02, we conclude that (pi(t) — Pj(t)) —> 0 as t —> + 00 , for all i, j G Af 
if G is rooted. The proof is complete. □ 

Theorem [2] gives a solution to the synchronization problem of the class of nonlinear systems 
([8]) with relaxed communication requirements. In fact, each agent needs to send its information 
to its prescribed neighbors only at some instants of time. This information transfer is also subject 
to constraints inherent to the communication channels such as irregular communication delays 
and packet loss. An important feature of the above result is that it gives sufficient conditions for 
synchronization, given in (1181) . that are topology-independent and can be easily satisfied with an 
appropriate choice of the control gains. Notice that the constant h* := ( k*T + h) can be easily 
estimated in practice, and is simply defined as the maximum blackout interval of time an individual 
agent does not receive information from each one of its neighbors. Then, the control gains, namely 
h 1 ? and A*, can be freely selected to satisfy G9; in particular, the variable pi can be made arbitrarily 
large, which is advantageous in the case where the parameter h* is roughly or over estimated. On 
the other hand, condition (flSl) is equivalent to 0 < h* < tj(m* — 1), which specifies the maximal 
allowable time interval during which each agent can run its control algorithm without receiving new 
information from its neighbors. This allowable interval of time does not rely on some centralized 
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information on the interconnection topology between the systems and can be made arbitrarily large. 
It should be also pointed out that the results of Theorem [2] are obtained under mild assumptions 
on the interconnection graph Q. In this regard, note that condition (TTSl) is imposed for all agents 
i £ Af\ where, in view of the assumptions on Q, the set Af\ Af* contains at most one element. 

Remark 2. The term eijft) in (fThl) can be selected in different ways using the estimates of the 
desired velocity of the i-th agent. The choices other than (El include: 

eij{t) := f v di {[s/T\)ds, (31) 

dkij (t)T 


and 


f-ij(t) := / v di (s)ds, 

JkiAW 


(32) 


■= v di (kij(t)T) ■ (t - kij(t)T). 


(33) 


In view of Proposition |7J any of the choices (fT71) . (131D - (1551) can be used for our purposes. 

The control scheme in Theorem [2] can be applied to the case where the desired velocity is 
available to all systems, i.e., £ = A f. In this case, the observer CT- dTTl) is not needed and the 
following result, which can be shown following similar arguments as the proof of Theorem [2j is 
valid. 


Corollary 1. Consider the network of n-systems ©, where the interconnection topology is de¬ 
scribed by a directed graph Q and the communication process between the systems satisfies Assump¬ 
tion QJ Suppose that each system is controlled by a control law Ui satisfying Assumption @ where 
the corresponding reference velocity v ri (t) is defined in ©, where v di = v d , i £ Af, and gift) is 
obtained from (fl4ll - (fT5l) with p^ ft) '■= Pj(kij(t)T) + v d ■ (t — kij(t)T), ( j,i ) £ £, and hij ft) is given 
in 0. Let the control gains satisfy El- Then, Vift) —> v d and fpift) — Pjft)) —> 0 as t —t +oo for 
all i,j £ Af and for arbitrary initial conditions if Q contains a spanning tree. 


5 Application to Euler-Lagrange systems 

In this section, we apply the proposed approach to the class of fully-actuated heterogeneous Euler- 
Lagrange systems. The systems dynamics are given by 

Pi = Vi 

Vi = Mi (pi) fUi - Cifpi , Vi)vi - Gi(pi )), 

for i £ Af, where pi £ R m is the vector of generalized configuration coordinates, Ui £ R m is the 
vector of torques associated with the i th system, Mifpf), Cifpi, Vi)vi, and Gifpf) are the inertia 
matrix, the vector of centrifugal/coriolis forces, and the vector of potential forces, respectively. The 
inertia matrices Mi(pi) are symmetrical and positive definite uniformly with respect to pi. Other 
common properties of Euler-Lagrange systems m are as follows: 

P.l The matrix Mifpf) — 2 Cifpi, vf) is skew symmetric. 

P.2 There exists k Ci > 0 such that \Ci(pi,x)y\ < fc Ci |x|-|y| holds for all pi, x, y £ R m . In addition, 
Mifpi) and Gifpf) are bounded uniformly with respect to pi. 

P.3 Each system in (1341) admits a linear parametrization of the form Mi(pi)£i + Cifpi, Vi)Xi + 
Gifpi) = Yi{pi,Vi,Xi,Xi)di, where Yifpi,Vi,Xi,ii) is a known regressor matrix and Oi £ M. k is 
the constant vector of the system’s parameters. 
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We assume that the systems are subject to model uncertainties; the parameters 0* in P.3 are 
unknown. We aim to achieve Objectives Q] and [2] for the Euler-Lagrange systems (|34p under a 
directed interconnection graph Q and the communication constraints described in Section 13.21 
For this purpose, we consider the following control input in (1341) 

Ui — Y[ {jpi , Vi , V ri , V Ti ]0i k^ Cj, (35) 

= hi j Yj f ('Pi , Vi, V Ti , V Ti ) d, (311) 

for i G Af, where the matrix LL is symmetric positive definite, Y) is defined in P.3, &i G is an 
estimate of the parameters, k] > 0, and et = (vi — v ri ) with 


v n 



Vi 


= Vi+ Vdi , 

= -kiVi - Ai(pi - 4>i) 

= -V ’i + ^ ^ £"=1 (i) 

= rji = 0, iGA/'\A/' # , 


i G A/" 8 , 


(37) 

(38) 

(39) 


where p^(t) := pj{kij(t)T)+Vd } {kij{i))-{t — kij{t)T), the control gains are defined as in Theorem^ 
kij(t) is defined in Q. v,i. is obtained from (fl2l) - (fl3l) with the discrete-time observer 


Vdi(cr) = v d , i G £, (40) 

Vdi {cr + 1) = ■ * ' ■ ^(cr) i&J 7 , (41) 

I I j G Ni(a) 

where id^ (cr) is given in m and Ni(cr) is defined after m ■ Then, the following result is valid. 

Corollary 2. Consider the network of n Euler-Lagrange systems (13411 interconnected according to 
Q and suppose that Assumption]] holds. For each system, let the control input be given in © with 
©-CUD. and suppose condition © is satisfied. Then, Objective 1 and Objective 2 are achieved 
under the conditions on the interconnection graph Q given in Theorem [H 

The proof of this result follows from Theorem [2] by noting that v ri G C 1 and v n is well de¬ 
fined, and the control law (155l) - (156l) is the standard adaptive control scheme proposed in gD] 
for Euler-Lagrange systems that satisfies Assumption [2] In fact, using the Lyapunov function 
Vi = |(el" Mi(pi)ei + 6 i Ur 1 6i), with 6t = ( 0i — Of), one can show that et G £2 H £<*, and 6i G £ 00 , 
leading to the first point in Assumption [2] Also, properties P.2 and P.3 guarantee that e, G 
if v ri , v ri G £00 • Then, invoking Barbalat Lemma, one can conclude that if v ri , v Ti are uniformly 
bounded, then e, (t) —> 0 as t —> + 00 , which is the second point in Assumption [2] 

Remark 3. The control scheme in Corollary]] extends the relevant literature dealing with Euler- 
Lagrange systems with communication constraints f8j fT7i[T9ll20ll26])27l for instance] to the case where 
a non-zero final velocity is assigned to the team, the communication between agents is intermittent 
and subject to varying delays and possible packet loss, and under a directed interconnection graph 
that contains a spanning tree. In addition, this control scheme extends the work in m to the case 
of intermittent and delayed communication without using a centralized information on the intercon¬ 
nection topology. Note that in JSjj, Objective]]is achieved, in the case of delay-free continuous-time 
communication between agents, under some topology-dependent conditions. 
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6 Simulation Results 

We provide in this section simulation results for the example in Section [5] Specifically, we consider 
a network of ten Lagrangian systems; A f = {1,..., 10}, modeled by equations (1M1) with m = 2, 
Pi ■= ( Ph,Pi 2 ) T > v i '■= K,tO T . and 


f 0 1 + 202 cos (p i2 ) 0 3 + 02 cos (p i2 ) 

V 03 + 02 COS (p i2 ) 0 3 


/ —02 sm(pi 2 )vi 2 - 9 2 sm(p i2 )(v il +v i2 ) \ 
\ 02 sinfejw^ 0 ) 


r( \_( 005 COS^iJ + 00 4 COs(p» 1 +Pi 2 ) 

11 V gd 4 cos(p h +Pi 2 ) 

where g = 9.81 m/sec 2 , the variables 6k, k = 1,. .., 5, are given as: 0i = (mi2^ + m 2 (l\ + 2 2 2 ) + h + 
h), 6 2 = m 2 l\l c2 , 0 3 = (m 2 2 2 2 + I 2 ), 04 = ni 2 l C2 , and 0 5 = (mil c i + m 2 h), with mi = m 2 = 1 kg, 
h = = 0.5 m, l c 1 = 2 C 2 = 0.25 m, and I\ = I 2 = 0.1 kg/rn 2 . The parametrization satisfying 

property P.3 for each system is given as: Yi(jpi, v,, x,, x,) = [Yi jk \ £ K 2x5 , with Yi 11 = x^, 1} 12 = 
cos(pi 2 )(2i;j 1 + x i2 ) - sm(pi 2 )(x il v i2 + x^^ + v i2 )), Y il3 = x ia , Y hi = Y i2i = gcos{ l pi + p i2 ), 
Y il5 = gcos( l pi), Y i21 = Y i25 = 0, F i2 2 = x it cos(p l2 ) + x il v i3 sin(p i2 ), Y i23 = x ix + x i2 , for any 
Xi := (a;i 1 ,Xi 2 ) T . The vector of estimated parameters is 0 = (0i, 02, 03 , 04 , 0s) T - 

The systems in the network are interconnected according to the directed graph Q given in Fig. 
[IJ with the node labeled 4 being one of its roots; r = 4. For the communication process, we set the 
sampling period T = 0.1 sec, which means that for each (j, i) £ £ 1 agent j can send its information 
to agent i only at instants /cT, k £ Z+. The delays and packet dropout are generated for each 
communication link as follows. For each (j, i) £ £, and at each instant kT , we pick the information 
p j(k), where k £ Z is randomly selected in the interval [k — 10, k\. This information is then delayed 
by r £ [0.15,0.25] sec and considered as received by agent i. Due to the random choice of k, a 
simple logic is implemented to avoid sending information at a future k with the same k. This way, 
the parameter h* is estimated to be 1.3 sec, the variable kij(t), for all t > 0, and the set A/}(cr) can 
be easily obtained, the information of agent j at some instants kT are lost (not submitted), and the 
information received by agent i is randomly delayed. The intermittent nature of the communication 
process as well as varying communication delays and packet dropout are illustrated in Fig. (2j which 
shows the received discrete-time signal <p(k) when the signal </>(£) = 2 sin(t) is sent according to the 
communication process described above. 

We implement the control scheme developed for Euler-Lagrange systems in Section [5] First, 
we consider the case where C = {1,4}, which indicates that the systems labeled 1 and 4 are the 





Figure 1: interconnection graph Q. Figure 2: Input and output of a communication chan- 

nel, (j){t ) = 2sin(t). 
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only systems having access to the desired velocity given by Vd = (0.3,0.6) T rad/sec. The observer 
(l9l)- (fTT1l is updated at T, and the control gains are selected as: k 1 ’ = kf = 2, IT = 0.31s, kf = 10, 
Aj = 13, kf = Note that this choice of the gains satisfies condition (fI51) with pt = y/Xf. The 

weights of the communication links of Q w , which is the same as Q with assigned weights on its links, 
are set such that k,; = 1. 

Fig.[3]and Fig.[]]illustrate the relative positions and relative velocities defined as pu = {pu^ ,pu 2 ) T 
(Pi ~ Pi)> for i = 2,..., n, and v\j = , r'ij 2 ) T := (Ti — Vj) for j = d, 2,..., n, where subscript 

l d’ is used for the desired velocity. It is clear that all agents synchronize their positions and veloc¬ 
ities with the desired velocity. The output of the discrete-time observer is given in Fig. [5l with 
Vdi = (vd i[L , i>d i2 ) T , where it can be seen that the desired velocity estimate of each agent converges 
to the desired velocity available to the leader agents. 




time (sec) 


Figure 3: Relative position vectors in the case of L = {1,4}. 




Figure 4: Relative velocity vectors in the case of C = {1,4}. 

Next, we consider the case where C = 0. Using the same above control parameters, the obtained 
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Figure 5: Desired velocity estimates in the case of C, = {1,4}. 

results are shown in Fig. [6][8] where v\j is defined for j = 2,..., n. These figures show that all systems 
synchronize their positions, and their velocities converge to the final velocity dictated by the output 
of the discrete-time velocity estimator. 




Figure 6: Relative position vectors in the case of C = 0. 


7 Conclusion 

We addressed the synchronization problem of second-order nonlinear multi-agent systems intercon¬ 
nected under directed graphs. Using the small-gain framework, we proposed a distributed control 
algorithm that achieves position synchronization in the presence of communication constraints. In 
contrast to the available relevant literature, the proposed approach guarantees that all agents veloc¬ 
ities match a desired velocity available to only some leaders (or a final velocity agreed upon by all 
agents in the leaderless case), while information exchange between neighboring agents is allowed at 
irregular discrete time-intervals in the presence of irregular time-delays and possible packet loss. In 
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Figure 7: Relative velocity vectors in the case of £ = 0. 


fact, we proved that synchronization is still achieved even if each agent in the team runs its control 
algorithm without receiving any information from its neighbors during some allowable intervals of 
time. The conditions for synchronization derived in this paper can be satisfied by an appropriate 
choice of the control gains. Future research will consider the extension of this work to the case of 
variable desired velocity. 


A Proof of Proposition [T] 

Consider the consensus algorithm ©-dD. The interaction between the agents in the system ©- 
ED is described by a directed graph Q s = (A f, £ s ), which can formally be obtained from the graph 
Q by modifying some of the its links, as follows: (1) removing the incoming arcs to each leader node 
(or agent), (2) adding a directed link from any leader node to any other leader node, and (3) adding 
a self arc to each node in the graph. It is straightforward to verify that, if the directed graph Q 
is rooted at r £ £, then Q s is also rooted at r £ C. In the case of no leaders {C = 0), the above 
modifications reduce to adding a self arc to each node; in this case, it is trivial that Q s is rooted if 
Q is rooted. 

In view of the above discussion, the consensus algorithm ®-ED can be formally written as 


v di (cr + 1) = 


ran ? 

1 ' 1 jeN^a) 


O)), 


for all i £ A/”, where 


Ni(a) = 


Ni{a) 

£ 


for ie J, 
for i £ £, 


(42) 


(43) 


and f^’^(tj) is a delay that takes some integer value at aT and, in view of Assumption Q] and ©, 
satisfies r^(a) < h% := \h*/T] for all a = 0,1,.. .. Note that f^\a) = 0, and f^V) = 0 for 
all a if i,j £ C. 

Let Gs{cr) = (N,£ s (cr)), where ( j,i ) £ £ s (cr) only if j £ Ni(a), which defines the set of those 
agents whose information is used in the update rule of agent i at instants aT. It is clear that G s (&) 
is a directed graph with at most one directed link connecting each ordered pair of distinct nodes 
and with exactly one self arc at each node. 
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Figure 8: Desired velocity estimates in the case of £ = 0. 


According to Theorem 2 in [42], the states of (14211 satisfy Vdi(cr) —> v c exponentially as a —» +oo, 
i £ J\f, for some v c £ R m , if the sequence of graphs Q s ( 0), (7 S (1),... is repeatedly jointly rooted. 
We claim that the latter condition is satisfied under the assumptions of Proposition [Q To show 
this, pick an arbitrary p £ 7L+ and consider the composition of graphs Gs(p) = Gs{p + h* a ) o Q s {p + 
ft* — 1) o ... o G s {p). Since G s {p) contains self arcs on each node, for all p, the edges of Gs{p + ft*), 
G s {p + ft* - 1), ..., Qs{p ) are also edges in Q s {p). 

Consider first the case where C ^ 0 and G is rooted at r £ C. Equation (l43l) implies that G s {p) 
contains a directed link from any leader node to any other leader node, for all p. In addition, in view 
of the definition of Gs = (J\f, £ s ), Assumption [T| implies that for each i £ T and (j, i) £ £ s , j ^ i, 
the information v^. is successfully delivered to agent i at least once per ft* := |~ft*/T] sampling 
periods. Therefore, it can be verified that, for each i £ T . if (j, i) is an edge in Gs, then (j, i) is also 
an edge in the composition of graphs Gs{p)- In fact, if (j. i) £ £ s , the definition of Ni(cr) implies 
that (j, i) is an edge in at least one of the graphs G s {p), G s {p + 1), • • G s (p + h%). Consequently, 
one can conclude that all the edges of Gs are also edges in the composition of graphs Gs{p), and 
therefore, Gs{p) is rooted at r £ C since Gs is rooted at r £ C. As a result, the sequence of graphs 
Gs (0), Gs (1), is repeatedly jointly rooted. Similar arguments can be used to show that Gs{p) is 
rooted in the case where C = 0 and G contains a spanning tree, and hence rooted. 

Now, since the states of the leaders are fixed, i.e., Ai, (c) = Vd for all a and i £ £, one can 
conclude that v c = v,j in the case where C ^ 0. The rest of the proof follows in view of the 
dynamics dsn, for i £ J 7 , which can be rewritten as A = — kfii — k^ii + fcf(Azi(LV^J) ~ v c), 
ti := Vdi — v c , and describe the dynamics of an asymptotically stable system with an exponentially 
convergent perturbation term. 
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